Easy Proofs of Some Borwein Algorithms for tt 



Jesus Guillera 



1 The Borweins' algorithms for n 

In 1987 Jonathan and Peter Borwein, inspired by the works of Ramanujan, derived 
many efficient algorithms for computing ir (see [3] and [4]). However their proofs 
are difficult for the nonspecialist in the theory of elliptic modular functions. We will 
see that by using only a formula of Gauss's and elementary algebra we are able to 
prove the correctness of two of them. The two algorithms we will consider are the 
quadratic algorithm: 

d = 1/V2, r = 1/2, 



1- Jl-d* 
1 + a/1 - dl 



r n+i — (1 + d n+ i) 2 r n — 2 n+1 d n+ i 



in which r n tends to 1/tt and in each iteration the number of exact digits is doubled, 
and the quartic algorithm: 



so 



l/</2, t = 1/2, 



S 4 



n 



tn+l = (1 + s n+l) 4 ^n — 2 2n+2 S n+ i (1 + S n +i + S^ +] /. 

in which t n tends to and the number of correct digits is quadrupled in each 
iteration. 
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2 A Gauss formula 

At the age of 14, C. F. Gauss, starting with two numbers a = a and b = b, 
considered the iteration 

a-n+i = — - — j b n+ i = \J a n b n , (1) 
and defined the arithmetic-geometric mean of a and b as the common limit 

M(a,b) = \ima n = \imb n . (2) 
In the year 1799, at the age of 22, by means of the elliptic integrals 

d9 r/ 2 cos 2 6d6 



U 2 n cos 2 9 + b 2 n sm 2 9 J ° a 2 n cos 2 9 + b 2 n sin 2 9 

he discovered the following wonderful formula pQ, pp. 87-102]: 

oo OM 2 

E^-^) = i--, (3) 

in which ao = 1, &o = 1/v2j a n and b n are given by (OQ), and M stands for the 
arithmetic-geometric mean of the numbers 1 and 1/Vz. Gauss's formula ([3]) was 
rediscovered in 1976 independently by E. Salamin and R. Brent (see [5] and [7]). 

3 Proof of the quadratic algorithm 

Our proof is based on the sequence 

u n ZU n k=n 

where ao = 1, &o = l/v2, a n and 6 n are defined as in (CD), and M = M(l,l/ v / 2), 
that is, M = lima n = limfe n . It is easy to see that limr n — 1/tt. Using relation ([3]) 
we also see that r = 1/2. In addition we can easily obtain 

r n a 2 n - r n+1 a 2 n+1 = 2 n " 1 (a 2 - 6 2 ), 

which leads to the recurrence 



To arrive at the Borweins' quadratic algorithm we write system (CD) in the form 



K _ 2 / a n 
a n \a n+ i 

and deduce from it the relation 



b n +i _ a n bn 

a n+l a n+l V a n 



(6) 



a n+l 



\ 



h 2 



b n+l 



(7) 



Finally we can see that if we define 



d n 



\ 



b 2 

1 - — 

at 



then 



and 



d n +i 



- 1 



a n -b n l-Jl-dl 



Wl On + ^n 1 + Jl - dl 



r n+1 = (1 + d n+ i) 2 r n - 2 n+1 rf n+1 . 
Substituting the values of ao and &o in (EI), we get the initial value do = l/y/2. 

4 Proof of the quartic algorithm 



(9) 
(10) 



From the relation ( ITOi) we have 

f 2n+2 = (1 + d2n+2) 2 r 2n+l ~ 2 2n+2 (i2n+2 

and 



T 2n+l — (1 + d2n+l)' 2 r 2n ~ 2 2n+1 (i2n+l • 



From © we get 



\ 



1 



l-d n+1 \ 2y/d~. 



n+l 



If we define s n = y/d^ then d 



1 + <i n+1 J 1 + <i n+ i 
2n+2 = ■s 2 +1 and (TTHT) implies that 



2n+l 



1 + S n+1 



'IF 



(12) 



(13) 



(14) 
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This allows us to rewrite ( fill) and ( 1121) as 

r 2n+2 = (1 + s 2 n+1 ) 2 r 2n+1 - 2 2n+2 s 2 n+1 (15) 

and 

r 2 „+, = (l + T ff-) 2 r 2 „ - 2« 2 T i^-. (16) 



Combining (1151) and (1161) and defining t n = r^n-, we obtain 

t n+1 = (1 + s n+ i) 4 t n - 2 2n+2 s n+ i(l + s n+ i + s 2 n+1 ). (17) 
On the other hand, from ([9]) we have 



i-^/i-4 



d2n+l — ; , (18) 

j l-si 



which together with (fT4"l) gives 

V " 1 + Sn + 1 

Finally we obtain 



1 - \ l - 4 

s n+1 = (19) 

i + di - st 



n 



5 Other algorithms for n 

Jonathan and Peter Borwein have proved the following cubic analog [2] of Gauss's 
formula (j3J): 

oo i jif2 

E^-«U= -- , (20) 
fc=0 13 7r 



where ao = 1, &o = 1/2 • y 18 — 6\/3. 



2 

and M = lima n = limfe„. 



fln+l — ~ , n+ i — \ , {Zi) 
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Rewriting system ( I2T1) in a way that is analogous to that used in ([6]) and defining 
e n+ i = — , r n = ~ + £ 3 fc (a^ - a 2 fe+1 ), (22) 

we recover the following Borwein cubic algorithm for 7r: 

e = \/3, r = 1/3, 

_ 3 

1 + y/8 - (e n - l) 3 

r n+l — e n+l r n ~ 3 (e n+1 — 1), 

in which r n tends cubically to 1/tt. In [6] we also find the following quartic analog 
of Gauss's formula (J3J): 

£/( a i " <i) = ^ IF' (23) 



where do = 1, &o = V 12\/2 — 16 



_ a w + fe n , 4 j a n b n (a 2 n + bl) 
a-n+i — — 2 — ' n+1 ~~ V 2 ' 

and M = lima n = lim6 n . 

In the same way as before, using this quartic analog and defining 

a 1 A/f 4 4 00 

en + i = ^, r n =^ + ^£4*(a -a* +1 ), (25) 
a„+i vr a 4 3a* fc=n 

we can obtain the following algorithm for it: 

e = a/2, r = 1/3, 
2 

1 + dl - (e n - l) 4 



, 4 n+1 



?"n+i — e n+1 r n ( e n+l 1)) 



3 

in which r n tends quartically to 
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